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Abstract. The Lorentz Integral Transform method (LIT) is extended to pion photoproduction in the A- 
resonance region. The main focus lies on the solution of the conceptual difficulties which arise if energy 
dependent operators for nucleon resonance excitations are considered. In order to demonstrate the appli- 
cability of our approach, we calculate the inclusive cross section for 7r°-photoproduction off the deuteron 
within a simple pure resonance model. 
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1 Introduction 

The Lorentz Integral Transform method (LIT) pQ has 
been proven to be a powerful technique for calculating in- 
clusive (see e.g. OOUTJ] and references therein) as well as 
exclusive |^]|^][7) photoreaction cross sections with com- 
plete inclusion of final state interaction (FSI) without cal- 
culating final continuum wave functions. Recently, the tech- 
nique has also been extended to electroweak processes [8]. 

This success has motivated the extension of the LIT 
method to pionproduction processes on light nuclei A > 3, 
. where existing approaches piTU] call for considerable im- 
provements. In the case of 3 He, for example, a conven- 
tional treatment of FSI would imply a four-body Faddeev- 
Yakubovsky treatment JI] of the final state which is very 
complicated. With respect to present experimental pro- 
grams to study electromagnetic meson production on light 
nuclei, e.g. at MAMI in Mainz, more sophisticated calcu- 
lations of such reactions are certainly needed in the near 
future. 

Recently, the LIT has been applied to inclusive pion 
photoproduction on the deuteron as the simplest possible 
nuclear target As a first step, only the near thresh- 
old region has been considered, where only the dominant 
Kroll-Ruderman |13| as production operator was included, 
and results comparable to traditional approaches were a- 
chicvcd. 

In the present work we want to extend this approach 
to higher photon energies into the Zl(1232)-resonance re- 
gion. As we will see below, one cannot proceed naively 
in a straightforward manner, because of conceptual prob- 
lems which arise from the energy dependence associated 
with the resonance contribution to the elementary pro- 



duction operator. We will first give in section [21 a brief 
outline of the LIT approach for energy independent tran- 
sition operators. The problem of the standard LIT method 
for electromagnetic particle production via resonances is 
discussed in section |21 where we also present a formal so- 
lution. As a test case, we consider in section 0] a simple 
model for 7r°-production on the deuteron in the Z\-region. 
The corresponding results, together with a summary and 
an outlook, are presented in section |3J 



2 The Lorentz Integral Transform method 
for energy independent transition operators 

In this section we review briefly the LIT method for inclu- 
sive reactions. The central quantity is the response func- 
tion 

i?H = J d!P> \(¥f\0\¥)\ 3 5(E f -Eo-u), (1) 

where O is an operator describing the transition from the 
ground state with energy Eq, to final states <Py, with 
energies -E 1 /, in the specific process under consideration. In 
the LIT approach, the response function is not calculated 
directly. Rather, one first introduces an integral transform 
of the response function R(u>) by 

L(a)= Tdu R{UJ \ 2 , (2) 

Ju th {Eo + W - <7R) Z + CTf 

where a = ctr, + i o\ with <j\ ^ 0. Furthermore, uj t h denotes 
the reaction threshold. 
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Inserting the response function from eq. l|T]l. one can 
rewrite eq. as follows 

L(a) = J dWd&f 5{E S - W) 

m — - — Wf)(&f\ — - — o m , 

= JdV f (V\Ol^—;\V f )(V f \^-^0\V). (3) 

Now the Schrodinger equation H\^f) = Ef\iPf) can be 
used to replace Ef with the hamiltonian H of the given 
system 

L(a) = / m f (*\ j^—t^fWA-E^—: O \W) . (4) 



<H-a 



In this step it is essential that the operator O is energy 
independent. By using the completeness relation 

^ d# f \V f )(V f \ = t, (5) 

one obtains finally 

L(a) = (if | O t (H - a*)- 1 (H - a)- 1 O \¥) 

= ($\V), (6) 

where one has introduced the so-called Lorentz state 

\$) = (H-a)- 1 0\V), (7) 

which obeys an inhomogcncous differential equation 



(H-a)\&)= 0\<F) 



(8) 



and which is bound at infinity. Recalling the fact that H 
is hermitean and therefore has only real eigenvalues, this 
feature guarantees a unique solution of (JSJ because the 
corresponding homogeneous equation has only the trivial 
solution. Since the source at the right hand side of (jHJ is 
localized and Imfer} ^ 0, the asymptotic behaviour of & 
at infinity is bound-state like. Thus the evaluation of L{a) 
avoids the explicit calculation of the continuum states <f/ 
but still includes the complete final state interaction. In 
the final step, the desired response function R is obtained 
from L(a) by an appropriate inversion method, see |14j 
concerning further details. 



3 The Lorentz Integral Transform method for 
meson production processes 

In the foregoing derivation an essential assumption was 
that the transition operator O is energy independent. This 
is, for example, the case for dipole absorption in the long 
wave length limit, but it is certainly no longer fulfilled for 
a retarded dipole operator. The same is true for the Kroll- 
Ruderman term in |12|. However, in both cases the energy 
dependence is smooth and weak. In principle, one could be 



tempted to replace then in the operator the energy by the 
corresponding hamiltonian. But then the equation for the 
Lorentz state becomes quite complicated and non-linear 
in the hamiltonian. In order to avoid this complication, 
another method has been devised [T2\ by treating the en- 
ergy in the operator as a parameter, fixed to some value 
e. One then determines a Lorentz transform L(a, e) as a 
function of this additional variable e. The inversion then 
yields a response function R(uj, e) from which the desired 
response function is obtained by setting 



R(u>) = R(u,uj). 



(9) 



However, this method fails for a resonance-like energy 
dependence which usually appears when certain degrees 
of freedom or states are projected out in favor of effective 
operators. This is the case in pion photoproduction where 
any realistic production operator contains resonance or 
pole contributions (see for example [L31ITC| ) describing in- 
termediate excitation of nucleon isobars. The problem, one 
encounters, can be illustrated already by the contribution 
of the nucleon-pole diagram to the production process 
on the two nucleon system depicted in Fig. which is 
an essential ingredient in pion production. Projecting out 
the intermediate iViV-system in time ordered perturba- 
tion theory, the corresponding amplitude has the following 
structure 

f 



0(E f ) 



(10) 



Ef + ie-Ha' 

where Hq denotes the free A^A^-hamiltonian. Recalling 
now the different steps in equation (TJJ) and JBJl, and apply- 
ing the LIT approach naively would mean to replace the 
energy Ef appearing in the propagator by the full hamil- 
tonian H of the final 7riViV-system. This, however, does 
not make any sense, because Ho and H act on different 
particle systems. Also the method of treating the energy 
as a parameter does not work, because in this case the 
effective operator becomes singular and thus the Lorentz 
state is not normalizable any more. 

The physical reason for the energy dependence of the 
contribution IjlOjl is quite obvious because the interme- 
diate JV TV-configuration has been projected out in favor 
of an effective operator. The same situation appears for 
the Z\-resonance contribution in Fig. \I\ which is also de- 
scribed by an effective, energy dependent operator. A way 
out of this dilemma is to enlarge the configuration space 
by including the previously projected-out states and thus 
avoiding the energy dependence of effective operators. 

In the present case, where a meson is produced, we 
have to include configurations with different numbers of 
particles, i.e. we have to construct a genuine Fock space 
description. We will first consider a general Fock space T 
consisting of N orthogonal subspaces labeled by I. The 
projectors Pi onto these subspaces fulfil 



N 



1 = ^2 Pi with PiP m = Si m Pi 



(11) 



1=1 



The full Fock space hamiltonian H = Hq + V consists 
of a kinetic part Hq which is diagonal and an interaction 
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Vim = PiVP m , with I, m S {1, N}, allowing transitions 
between the various subspaces. A transition operator, de- 
scribing an electromagnetic process, is generally expressed 
by a hermitean current operator in J 7 , which does not 
contain any resonance-like energy dependence, 



0= J2 e ^lrn 



With ft m = P l fP m , 



(12) 



where denotes the photon polarization vector. A re- 
maining trivial energy dependence via current structures 
depending on the photon momentum can be handled by 
the parameter method. 

In this Fock space, the LIT equation reads now 



[H - <j)\V) = 0\W) , Im{cr} ^ . 



(13) 



To solve this equation we use the method of resolvents in 
analogy to standard scattering theory, i.e. we introduce 



00 M = 



1 



tr-Ha 



G(a) = 



1 



a - H 



(14) 



Defining a conventional T-matrix by TGo = VQ, one then 
obtains the following standard equations for the Q- and 
T-matrices 



G = Go + GoTGo = Go + GoVG . 
T = V + VG T . 



(15) 
(16) 



A formal solution of the LIT-equation can be written as 

W) = -go\w) . (17) 

Since O is regular and Im{cr} ^ , the Lorentz state \fr is 
normalizable like a bound-state. 

4 Application to 7r°-photoproduction via the 
excitation of the Z\-resonance 

As a test case, we will consider now 7r°-photoproduction 
off the deuteron in the Z\-region. Then T comprises be- 
sides ir°NN- also NN- and iVZ\-states, so that nucleon- 
pole as well as Z\-pole contributions to the photoproduc- 
tion process can be taken into account. Using a compact 
matrix notation with respect to the three projectors P„, 
P N and P A for the ir°NN, NN- and AZ\-subspaces, re- 
spectively, we can cast kinetic and interaction parts of the 
hamiltonian H = T + V into the following form 



VnA Vntt 
VaA Va-k 
V-rrA V n7r 



(18) 



Transitions between different subspaces are described by 
the non-diagonal matrix elements of V . For example, V^n 
contains the 7r JV-vertex allowing transitions between the 
NN- and 7r°jV./V-states. Since we only want to demon- 
strate the applicability of our approach and do not aim 



at a quantitative description of the data, we restrict the 
interaction solely to Vak and Vnn, i-s. we use 



'Vnn 
V 10 Va-k 
V nA 



(19) 



The parametrization of the AN 7r°-vertex Va-k is taken 
from |17j . As electromagnetic current we consider here 
solely the dominant Ml-iVZ\-current j NA 



jan\ k ) 



2M 



N 



i (tan x k 



(20) 



with G°J[ N = 4.22. With these building blocks we can 
describe the dominant Z\-pole contribution to photopi- 
onproduction. With the approximations (|19|l and l|2()|) , the 
iVW-interaction Vnn contributes solely to the deuteron 
ground-state. For reasons of simplicity, we have used a 
pure S-w&ve Yamaguchi potential [18| with modern pa- 
rameters from PHI ■ For a shorter notation we label opera- 
tors A connecting the same subspace with only one index, 
i.e. Au = Ai. The Lorentz state can then be written as 
(see Fig. 

- W) = GO A nW) = (1a + Go-kV v a)GaOanW) -(21) 
The resolvent Ga fulfils the equation (see Fig. |3J 

Ga = GoA + GqaVa-kGq-kV-kaGa , (22) 

which allows one to rewrite the LIT into the following 
form 

L = -—1^{{^\O na GaOanW)}. (23) 

Va consists of a loop diagram (disconnected) and a gen- 
uine two-body one-pion exchange potential (connected), 
see Fig. [3] If the latter is neglected, we will refer to it as 
impulse approximation (IA) and denote the correspond- 
ing resolvent from here on as G l A - Therefore in impulse 
approximation L 1A is given by i|23|) where Ga is replaced 

by 

G l £ = T 1 T i r , (24) 

O — Ta - ^A\<T) 

with Sa as self-energy of the A 

Ea (p) = Va-k Go-k (f) V„a I disconnected ■ 

(25) 



5 Results 

The result for 7r°-photoproduction on the deuteron is shown 
in Fig. 21 Besides the IA according to (|24H with a vari- 
able cr-dependent Z\-mass, we have calculated in addition 
another IA with a static Z\-mass (M A = 1232 MeV) for 
a better comparison with the IA of [201 . Considering the 
small model differences to [2U] the agreement is quite satis- 
factory. While the IA is calculated without a partial wave 
decomposition, it is needed for the 7VZ\-interaction con- 
tribution generated by the one-pion exchange diagram in 
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Fig. El In order to fulfil convergence, we took into account 
all channels with a total angular momentum J < 5 of the 
7VZ\-system. One readily notes that the iVZ\-interaction 
has a sizeable influence near threshold leading to an en- 
hancement and is still moderate near the maximum result- 
ing in a slight upward shift of the position and lowering of 
the absolute size by about 8 %. For a realistic description 
of the cross section in the Z\-region one has to include 
besides non-resonant photoproduction contributions the 
neglected FSI, i.e. NN- and irN- interactions \'21\ . 

In conclusion, we have demonstrated how to extend 
the LIT method when energy dependent effective opera- 
tors are involved. It turns out that in such a case those 
states, which have been projected out in favor of effec- 
tive operators, have to be included explicitly in an ex- 
panded Fock space. With respect to the specific process 
considered here as an example, pion photoproduction on 
the deuteron, it is obvious that in the future the present 
model has to be extended towards a more realistic descrip- 
tion of this reaction. First of all, additional FSI as well 
as nonresonant Born contributions have to be included. 
Furthermore, since this method really pays-off for more 
complex systems with more than two nucleons, one has to 
consider meson-photoproduction on other light nuclei like 
3 He or 4 He. 

This work was supported by the Deutsche Forschungs- 
gemeinschaft (SFB 443). 
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Fig. 1. Nuclcon- and Z\-pole contributions to pion photopro- 
duction on the deuteron. 
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Fig. 2. Diagrammatic representation of $ in eq. I|21|l . 



References 



Qa = 



+ 



8. 
9. 
10. 

11. 

12. 

13. 
14. 



15. 



16. 



17. 



VD. Efros, W. Leidemann, and G. Orlandini, Phys. Lett. 
B338 130 (1994). 

S. Martinelli, H. Kamada, G. Orlandini, and W. Glockle, 
Phys. Rev. C 52 1778 (1995). 

S. Bacca, M. Marchisio, N. Barnea, W. Leidemann, and 

G. Orlandini, Phys. Rev Lett. 89, 052502 (2002). 

S. Bacca, H. Arenhovel, N. Barnea, W. Leidemann, and 

G. Orlandini, Phys. Lett. B 603 159 (2004). 
VD. Efros, Sov. J. Nucl. Phys. 41 949 (1985). 

A. La Piana, and W. Leidemann, Nucl. Phys. A677 423 
(2000). 

S. Quaglioni, W. Leidemann, G. Orlandini, N. Barnea, and 
VD. Efros, Phys. Rev. C 69, 044002 (2004). 
D. Gazit and N. Barnea, Phys. Rev. C 70, 048801 (2004). 
J. M. Laget, Phys. Rep. 69, 1 (1981) 

S. S. Kamalov, L. Tiator, and C. Bennhold, Phys. Rev. 
Lett. 75 1288 (1995). 

W. Glockle, The Quantum Mechanical Few-Body Prob- 
lem, Springer (1983). 

C. ReiB, W. Leidemann, G. Orlandini, and E.L. Tomusiak, 
Eur. Phys. J. A 17, 589 (2003). 

N.M. Kroll and M.A. Ruderman, Phys. Rev. 93 233 (1954). 

D. Andreasi, W. Leidemann, C. ReiB, and M. Schwamb, ac- 
cepted for publication in Eur. Phys. J. A (2005); preprint 
nucl-th/0503033 

I. Blomqvist and J.M. Laget, Nucl. Phys. A280 405 
(1977). 

D. Drechsel, O. Hanstein, S. S. Kamalov, and L. Tiator, 
Nucl. Phys. A645 145 (1999). 

H. Popping, P. U. Sauer, and X.-Z. Zang, Nucl. Phys. 
A474, 557 (1987). 



+ 



Qa + 



Qa 



Fig. 3. Diagrammatic representation of Qa in eq. 122H . 
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Fig. 4. Total cross section of 7r°-photoproduction on the 
deuteron. Notation of the curves: dashed: IA with constant 
Z\-mass; dashed-dotted: IA using eq. 112H : full: inclusion of 
N A- interaction; dotted: IA according to j20| (only resonance). 



